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Abstract 



. ^ The IGEX experiment has been operating enriched germanium detec- 

^ ■ tors in the Canfranc Underground Laboratory (Spain) in a search for the 

neutrinoless double decay of ^^Ge. The implementation of Pulse Shape 
Discrimination techniques to reduce the radioactive background is de- 
scribed in detail. This analysis has been applied to a fraction of the 
IGEX data, leading to a rejection of ~60 % of their background, in the 
region of interest (from 2 to 2.5 MeV), down to ^ 0.09 c/(keV kg y). 



1 Introduction 

The nuclear Double Beta Decay (DBD) is an unique laboratory to investigate 
the nature and properties of the neutrino ^ |2I . The neutrinoless decay mode, 
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if it exists, would provide an unambiguous evidence of the Majorana nature of 
the neutrino, its non-zero mass, and the non-conservation of lepton number. 
After the definitive confirmation that neutrinos have indeed non-zero mass, as 
the solar and atmospheric neutrino oscillation results imply [310], the neutrino- 
less Double Beta Decay (DBD) has become a most relevant subject of research 
because it is a process able to provide, in a relatively short time, the neutrino 
mass scale and its hierarchy pattern. The current best sensitivity limits of the 
effective Majorana electron neutrino mass -derived from the neutrinoless half-life 
lower bound- stands around {niy) < 300 — 1300 meV 0101 (the dispersion range 
is due to the uncertainties in the evaluation of the nuclear matrix elements), 
whereas most of the neutrino mass models indicate that the mass region where 
the Majorana nature of the neutrino would be resolved lies two orders of mag- 
nitude below (~ a few meV). To achieve the sensitivity required for such new 
objectives, it will require a large number of double beta emitter nuclei, a very 
low background and a sharp energy resolution in the Q-value region, and/or 
methods to disentangle signal from noise. 

A typical example of this type of search is IGEX. The International Ger- 
manium Experiment (IGEX) has been running in the Canfranc Underground 
Laboratory at a depth of 2450 m.w.e. in a search for the neutrinoless Double 
Beta Decay. Details and results of the experiment can be found in ref. [5]. 
Three Germanium detectors (RGl, RG2 and RG3), of ~2 kg each, enriched to 
86% in '^^Ge were used. 

This paper presents the development and application of one technique used 
to reduce part of the radioactive background by discriminating it from the ex- 
pected signal by comparison of the shape of the pulses (PSD) of both types of 
events. As an example, the method has been applied to the data recorded by 
two Ge detectors of the IGEX Collaboration, which has produced one of the 
two best current sensitivity limits for the Majorana neutrino mass parameter 
00. 

The rationale for PSD is quite simple: in large intrinsic Ge detectors, the 
charge carriers take 300 - 500 ns to reach their respective electrodes. These 
drift times are long enough for the current pulses to be recorded at a sufficient 
sampling rate. The current pulse contributions from electrons and holes are 
displacement currents, and therefore dependent on their instantaneous velocities 
and locations. Accordingly, events occurring at a single site (/3/9-decay events for 
example) have associated current pulse characteristics which reflect the position 
in the crystal where the event occurred. More importantly, these single-site 
events (SSE) frequently have pulse shapes that differ significantly from those 
due to the background events that produce electron-hole pairs at several sites by 
multi-Compton-scattering process, for example (the so-called Multi-Site Events 
(MSE)). Consequently, pulse-shape analysis can be used to distinguish between 
these two types of energy depositions: DBD events belong to the SSE class of 
events and will deposit energy at a single site in the detector while most of the 
background events belong to the MSE class of events and will deposit energy at 
several sites. 

The IGEX detectors have modified preamplifier electronics to route and 
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record the current pulses at a very early stage of preamplification, thus produc- 
ing unique high-bandwidth pulse shape signals. Furthermore, to develop PSD 
techniques it would be highly desirable to obtain an earlier signal, even before 
it passes through the few unavoidable electronic components at the first stage 
of the detector preamplifier, resembling as much as possible the displacement 
current of the detector. This allows the development of algorithms that do not 
depend strongly on the preamplifier electronics in use. To this end, the transfer 
function of the preamplifier and associated front-end stage has been measured 
for each detector. This allows the reconstruction of the displacement current 
and an easy comparison to computed pulse shapes. 

The paper is organized as follows: in Section |21 the properties of the digi- 
tized experimental pulses are shown and the work performed to understand and 
reproduce their shapes is presented. In Section^ the PSD method is described. 
Finally, Section 01 displays the results of this analysis for the IGEX data. 

2 Experimental and calculated pulses 

The data acquisition system of the IGEX experiment is based on standard NIM 
electronics, each Ge detector having an independent electronic chain. Pream- 
plifiers were modified for pulse shape analysis and each preamplifier fast-pulse 
output is routed to a LeCroy 9362 digital oscilloscope (800 MHz analog band- 
with). The digitized pulse signal covers a total time of 1 ^s using 500 points; it 
is worth noting that the time resolution, of about 100 ns (as inferred from the 
width of the peaked features or the fall time) , limits the ability to resolve nearby 
features in the pulse such as lobes or discontinuities characteristic of a multiple- 
site interaction signal. Figurenshows the main features of the digitized pulses. 
These output pulses are taken at the very first stage of the amplification chain, 
but even so, there is an unavoidable instrumental distortion due to the pream- 
plifier. This has been studied to determine its transfer function h{t). To take 
into account this distortion in the pulse shape analysis, either the calculated 
pulse i{t) is folded with the transfer function. 



or the experimental signal unfolded. The latter allows the recovery of some 
information lost because of the instrumental distortion 7 . The transfer function 
of the preamplifiers, h(t), depicted in figure |21 has been directly measured as 
the response of the preamplifier for a narrow (5-like signal. Studies were also 
made following several methods including analog simulation of the preamplifier 
circuits and the analysis of the shapes of selected populations of experimental 
pulses. It was observed that the shape of the rise and the fall of the folded pulses 
is the same that the shape of the transfer function for those events in which the 
energy is released in some particular regions of the crystal; in particular, the 
left side of the transfer function can be deduced by studying pulses of events 
in the inner and lower part of the crystal, while the right side is derived from 




(1) 



3 



those produced in the outer and upper region. Measurements and estimates for 
the transfer function are found to be in quite good agreement. 

The pulse shapes of the output signals can be reproduced numerically. An 
energy deposition in a Ge crystal produces a proportional number of electron- 
hole pairs, which move towards the electrodes. The induced current i, and 
consequently the electric pulse taken from the detector, is the sum of the con- 
tributions due to each type of charge carriers: 

z(<) = Zc(0+jhW (2) 
The current is calculated as 0: 

«o(h) = -go(h)^wWo(hj(^) (3) 

where is the so-called weighting field"*. In principle, the drift velocity v of 
charge carriers is proportional to the electric field E: v(r) = fiE(r), where /x 
is the mobility in the material medium. However, at high values of the electric 
field, a saturation velocity is reached. The dependence of the velocity on the 
electric field has been obtained empirically 



-(E) = T^^^i^T^ (4) 



The values commonly used for the parameters involved in expression 01 (the 
mobility /i, the numerical parameter /3 and the electric flield for saturation £^sat) 
are summarized in table ^ 10 . 



Table 1: Values commonly considered for the parameters involved in the empir- 
ical dependence of the velocity of the charge carriers on the electric field. 





electrons 


holes 


E,,t(V/cm) 


275 


210.5 


^ (cm^/Vs) 


36000 


42000 


P 


1.32 


1.36 



The solution for the electric field (i?("f*) — ^0(T^)) in the crystal is derived 
from the Poisson equation, applied for the depleted and non-depleted regions: 

VV(T^) = -Pi^, depleted (5) 

V'^(t){T') = 0, non - depleted (6) 

4'i'r^int) - H^cxt) = Vo (7) 

*The weighting field is the field that would appear in the crystal if a unity voltage was 
applied to the sensor electrode (assuming no impurity in the crystal) 
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The electric field depends on the geometry of the crystal (true coaxial or 
closed-end), the supplied voltage Vq and on the residual space charge density 
p. The parameters rint and Text correspond to the locations of the internal and 
external electrodes of the detector. A solution (depicted in figure |3Jl has been 
obtained for the cases of depleted and non-fuUy depleted, closed-end crystals of 
the IGEX experiment by using an iterative method of calculation. For a non- 
fuUy depleted detector, the boundary between the depleted and non-depleted 
regions in the crystal is unknown a priori. An invalid solution is obtained 
when using the boundary condition of eq. {Tj), because a region with a field 
having inverse direction appears (which means that this region is not depleted). 
Therefore, the boundary condition is changed by using the deduced limit for the 
depleted region instead of rjnt ; then, a new solution is found for the electric field 
and, consequently, a new limit for the non-depleted region is obtained. This 
procedure is repeated until the solution does not change significantly between 
two successive iterations. 

Once the electric field is known, the pulse shapes can be calculated. Some 
examples for SSE are shown in figure ^ Different radial positions have been 
considered in each plot for several vertical coordinates; the effect of the instru- 
mental distortion is shown on the right plots. To reproduce MSE pulses, the 
shapes due to each individual interaction should be properly added weighted 
by the fractional amount of energy released in the crystal, obtained by Monte 
Carlo simulation. 



3 PSD Method 

The PSD method we have used consists in counting the number of lobes of 
the pulses and rejecting those events having more than two significant lobes or 
peaks. A SSE pulse is expected to have at most two lobes, one due to electrons 
and the other due to holes. Experimental pulses are first unfolded using the 
transfer function of the preamplifier. Then, to detect lobes a " mexican-hat" 
filter F of the proper width cr is applied to the pulse. In fact, this filter is the 
second derivative of a gaussian: 

nO = ^xexp(-|,) (8) 

and the filtered signal has a null mean value where there is no lobe in the original 
signal and a peak where a lobe is present. Therefore, it is straightforward to 
reject all the events having more than two lobes. Figure |H shows the results of 
applying this method to four different pulses. 

This robust method is nearly model-independent. Its effectiveness has been 
evaluated on calculated SSE. A test population of 2000 SSE pulses was generated 
for each detector. The locations were randomly chosen, uniformly distributed 
in the volume of the crystals. Each calculated pulse was folded with the proper 
preamplifier transfer function, then scaled to unit height, and a variable amount 
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of gaussian noise was added, to reproduce experimental pulses of different ener- 
gies and noise levels. Finally, the number of lobes of each pulse was obtained by 
applying " mexican-hat" filters derived from gaussians of different widths. After 
the analysis, the fraction of pulses with three o more lobes was retained. This 
is an estimation of the expected number of misidentified SSE pulses with this 
technique. The results are shown in figure El When applied to a mixed set of 
SSE and MSE pulses, to obtain the best discrimination results, a narrow filter 
should be preferred over a broad filter, but it would produce a large fraction 
of misidentifications. Conversely, a safer, broader filter would not find as many 
lobes as a narrower filter, and its discrimination power would be smaller. As 
a compromise, we choose the filter with a characteristic width of 60 ns, thus 
keeping the misidentification error for the calculated SSE pulses under 5% for 
all pulses in the whole noise range considered. Notice that those MSE in which 
the separation of the interactions in the (r,z) plane is too small or the amounts 
of energy deposited at each point are extremely unbalanced will not be rejected 
by this technique. 

This method of counting the number of lobes of the pulses has been applied 
also to a ^^Na calibration spectrum and to a set of data taken following a large 
intrusion of radon in the shielding. Fig. [7| illustrates the reduction in the case 
of ^^Na spectrum for detector RG2 and Fig. |S1 the case of radon for detector 
RG3. A comparison of the results of the method for the cases of background, 
^^Na calibration and radon is shown in Fig. |51 

4 Results 

The background reduction technique described above has been applied to the 
IGEX data whose pulse shape was recorded (those events whose pulse shape 
was not available are conservatively considered SSE). Table |5] summarizes the 
results (exposure, background levels in the region of interest from 2 to 2.5 MeV 
before and after the PSD, rejection factor) for each detector. Figure [TUI shows 
the spectra before and after application of the PSD technique to the data of 
detectors RG2 and RG3. This method results in an efficient rejection leading 
to a background level (in the best case) of 0.10 c/(keV kg y) in detector RG2, 
as can be seen in Tabled This value can be considered as the background limit 
achieved with the present PSD technique. The overall final background level of 
the set of detectors RG2 and RG3 together turns out to be of 0.10 c/(keV kg 
y), in the region of interest. 

The PSD analysis has been applied to only 52.51 mole y out of the total 
116.75 mole y (8.87 kg y) accumulated in the IGEX experiment. In Tabled 
the IGEX data correponding to 8.87 kg.y in ^^Ge in the region between 2020 
and 2060 keV, in 2-keV bins, are given, with and without application of PSD 
(see also ref. |2])- The obtained half-fife lower bounds are > 1-13 x lO^^j/ 
for the complete data set and of T^^^ > 1.57 x lO^^y for the complete data set 
with application of PSD to 52.51 mole y. Accordingly, the upper limits on the 
neutrino mass parameter are 0.38-1.55 eV for the first data set and 0.33-1.31 eV 
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Table 2: Results of applying the PSD (exposure, background levels b in the 
2-2.5 MeV region before and after the discrimination and rejection factors). 





exposure 


b before 


b after 


rejection 
factor 




kg y 


c/(keV kg y) 


c/(keV kg y) 


(%) 


RG2 


2.75 


0.27 


0.10 


62.19 


RG3 


1.90 


0.26 


0.11 


57.61 


total 


4.65 


0.26 


0.10 


60.36 



for the second data set pi. The uncertainties originate from the spread in the 
values of the calculated nuclear matrix elements. 

5 Conclusions 

A Pulse Shape Discrimination technique to reject the radioactive background in 
the region in which the double beta decay signal is expected has been developed 
and applied to the data collected in the IGEX experiment, searching for the 
neutrinoless double beta decay of ^^Ge. A satisfactory understanding of the 
pulse shapes has been achieved. The method described in this paper is based 
on the counting of the number of lobes of the pulses, using a proper filter. It has 
provided a rejection of ^ 60 % of the events in the region of interest, accepting 
the criterion that those events having more than two lobes cannot be due to a 
double beta decay. Accordingly, the improved background levels provided by 
the PSD technique have allowed the improvement of the limits for the half-life 
of ^^Ge and consequently, the effective electron neutrino mass bound. 
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Figure 1: Main features of the digitized experimental pulses. 
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RG1 V=4800 volt iteration-1 RG2 V=3800 volt iteration-2 RG3 V=3800 volt iteration-4 




Figure 3: Electric field for detectors RGl, RG2 and RG3. The operation voltage 
and the number of iterations necessary to derive the field (see text) are shown 
for each detector. 
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Figure 4: Calculated pulse shapes for SSE in detector RG2. Different radial 
positions have been considered in each plot for several vertical coordinates h (H 
is the height of the crystal). On the right, the instrumental distortion has been 
taken into account. 



12 



setup#10 2090-18 E=2043.5 keV L=2 SSE setup#10 3042-7 E=2022.6 keV U3 MSB 




setup#10 3080-360 E=2045.4 keV L=2 SSE setup#10 2037-98 E=2050.0 keV L=3 MSE 




time time 



Figure 5: Examples of the effect of applying the " mexican-hat" filter to detect 
significant lobes in the digitized pulses. Events on the left are accepted (having 
two lobes) while those on the right are rejected (having three lobes). 
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C haraderistic filter width (ns) C liaracteristic filter width (ns) 

Figure 6: Fraction of misidentified calculated SSE pulses as a function of the 
characteristic widht of the analyzing filter and the amount of noise added to the 
pulse (see text). 
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Figure 7: Energy spectra before and after the PSD based on the counting of the 
number of lobes for a ^^Na calibration and for detector RG2. A zoom of the 
plot at the top is shown at the bottom. 
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Figure 8: Energy spectra before and after the PSD based on the counting of 
the number of lobes for radon data and for detector RG3. A zoom of the plot 
at the top is shown at the bottom. 



16 




500 1000 1500 2000 2500 

E (keV) 



Figure 9: Comparison of the percentage of SSE events identified for a ^^Na 
calibration, radon data and background, for both detectors RG2 (top) and RG3 
(bottom). 
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Figure 10: Background spectra before and after the PSD based on the counting 
of the number of lobes for detectors RG2 (top) and RG3 (bottom). 
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